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Abstract. We have derived the dust emissivity in the 
Far-Infrared (FIR) using data available in the literature. 
We use two wavelength dependences derived from spectra 
of Galactic FIR emission (Reach et al. 1995). A value for 
the emissivity, normalised to the extinction efficiency in 
the V band, has been retrieved from maps of Galactic 
FIR emission, dust temperature and extinction (Schlegel 



et al. 1998) 



Our results are similar to other measurements in the 
Galaxy but only marginally consistent with the widely 
quoted values of Hildebrand ( 1983 ) derived on one reflec- 
tion nebula. The discrepancy with measurements on other 



reflection nebulae (Casey |l99l| ) is higher and suggests a 
different grain composition in these environments with re- 
spect to the diffuse interstellar medium. 

We measure dust masses for a sample of six spiral 
galaxies with FIR observations and obtain gas-to-dust ra- 
tios close to the Galactic value. 
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1. Introduction 

Assessing the quantity of dust in spiral galaxies is of 
primary importance in both understanding the intrinsic 
properties of galaxies themselves and interpreting obser- 
vations of the distant universe: large quantities of dust 
can modify the optical appearance of galactic structures 
like spiral arms (Trewhella 1998); if the distribution of 
dust is extended, a large fraction of the radiation from 
the distant universe can be blocked (Ostriker & Heisler 



1984); star formation as determined from UV fluxes could 



be severely underestimated thus altering our knowledge 
of the star formation history of the universe (Hughes et 
al. |199§). 

Dust mass can be retrieved from extinction or from 
emission in the FIR. In the former case information about 
the star-dust relative geometry is needed and the method 



can only be applied to nearby edge-on galaxies, where the 
dust distribution can be inferred from extinction features 



(Xilouris et al. 1997, 1998). In the latter case there are 
no such limitations, and the wealth of data in the FIR 
and Sub-mm from instruments like the Sub-mm camera 
SCUBA and from the satellites ISO and COBE, can be 
used to measure dust mass. Unfortunately, the determina- 
tion of dust mass is entangled with that of dust temper- 
ature and they both rely on knowledge of the dust emis- 



sivity (Hildebrand 1983), the form of which is currently 
highly uncertain. 

The emissivity (or emission efficiency, i.e. the ratio be- 
tween the emission cross section and the geometric cross 
section), <3 m(A) is usually described by a function of the 
form 



Qcm(A) — Qcm(Ao) 



(1) 



where Q e m(Ao) is the value of the emissivity at the refer- 
ence wavelength Ao, and f3 is the wavelength dependence 
index. 

While a value (3 = 1 seems to be plausible for A < 
100 (im (Hildebrand 1983; Rowan-Robinson 1992] ), there 
is observational evidence for a steeper emissivity at longer 
wavelengths. The difference in emissivity is not unex- 
pected, since emission in the Mid-Infrared (25-60 /urn) is 
dominated by transiently heated grains, whil e at A > 100 
grains emit at thermal equilibrium (Whittet 1992). Sub 
mm observations of spiral galaxies (Bianchi et al 



1998 



Alton et al. |1998b) show that it is not possible to use an 



emissivity with = 1 to fit the 450 and 850 /an emission. 



Reach et al. (1995) came to a similar conclusion. They 
used the spectrum of the Galactic plane observed by the 
spectrophotometer FIRAS on board the satellite COBE, 
to find that the data are well fitted by an emissivity: 



>(A)« 



\- 



l + (Ai/A) e 



1/6' 



(2) 
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for the range 100 /im to 1 cm. Eq. (0) behaves like (fil) 
with (3 = 1 at small A (A <C Ai) and (3 = 2 at large A 
(A > Ai) (they set A : = 200-^m). 
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Masi et al. ( 1995| ) measure a value (3 = 1.54 by fitting a 
single temperature grey-body spectrum to Galactic plane 
data in four bands between 0.5 and 2-mm taken by the 
balloon born telescope ARGO. Reach et al. (1995) sug- 
gest that a single temperature fit may bias towards lower 



values of (3 (see also Wright et al. 1991); over the whole 
FIRAS spectral range, a two temperature grey-body with 
(3 = 2 at large A provides a significantly better fit than a 
single temperature spectrum with f3 rs 1.5. At long wave- 
lengths theoretical calculations for crystalline subs tance s 
constrain f3 to be an even integer number (Wright 1993|) . 
For amorphous mater ials (3 depends on the temperature: 
Agladze et al. ( 1996 ) find 1.2 < (3 < 2 for amorphous 
silicate grains at a temperature of 20 K. 

A value for the emissivity at a specific wavelength 
Qom(Ao) normalised to the extinction efficiency in the op- 
tical can be determined by carrying out an energy bal- 
ance in a reflection nebula, comparing the energy absorbed 
from the central star with the FIR output from the sur- 
rounding dust. Alternatively, the extinction measured to- 
ward the star can be directly compared to the optical 
depth in the FIR (Whitcomb et al. [l98lt Hildebrand ~~ 



1983 



Casey |1991| ). These methods are complicated by the un- 
known nebular geometry and by temperature gradients in 
the dust; as an example, Casey (1991) found that the ex- 
tinction method usually retrieves higher values than the 
energy balance. 

In this paper we use the extinction method comparing 
the Galactic extinction to FIR emission: in this case the 
same column density of dust is responsible both for emis- 
sion and extinction and a reliable result can be obtained. 

2. The method 



Schlegel et al. (^998[ hereafter SFD) have presented a new 
map of Galactic extinction. After removing emission from 
zodiacal light and a cosmic infrared background, they have 
combined the 100 /im map of Galactic emission taken by 
the DIRBE experiment on board the COBE satellite with 
the 100 /im large-area ISSA map from satellite IRAS, to 
produce a map of Galactic emission with the quality cal- 
ibration of DIRBE and the high resolution of IRAS. The 
dust temperature has been retrieved using the DIRBE 
maps at 100 /im and 240 /im assuming /3=2. Knowing the 
temperature, the 100 /im map has been converted into a 
dust column density map and subsequently calibrated to 
E(B-V) using colours and Mg 2 -index of elliptical galaxies. 
We would like to stress that the colour excess has been 
derived from the 100 /im emission without any assump- 
tion about the value of the emissivity at any wavelength. 
Moreover, the choice of (3 does not affect significantly their 
results: when /3=1.5 is used, the dust column density map 
varies only of 1%, aside from an overall multiplicative fac- 
tor that is taken account of when calibrating with the 
colour excess. We have accessed the electronic distribu- 
tion of this remarkable dataset to retrieve the 9.5'/pixel 



maps of the intensity at 100 /im, 1(100 /im), the tempera- 
ture and the colour excess E(B-V) for the north and south 
Galactic hemispheres. 

When the same dust grains are responsible for emission 
and extinction, the ratio between the extinction coefficient 
in the V-band and the emissivity at 100 /im is equivalent 
to the ratio of the optical depths 



t(V) 



t(V) 



^100 /im) r(100/im)' 



(3) 



The above formula is correct if all of the dust grains are 
identical. In a mixture of grains of different sizes and ma- 
terials, the ratio of emissivities in Eq. (0) can still be re- 
garded as a mean value characteristic of diffuse galactic 
dust, if the dust composition is assumed to be the same 
on any line of sight. 

The optical depth at 100 /im, in the optically thin case, 
is measured using 



t(100 /im) = 



7(100 /im) 
B{100 /rni,T d y 



(4) 



where 5(100 /im, T^) is the value of the Planck function 
at 100/im for a dust temperature T^, both the intensity 
7(100 /im) and T& coming from the maps of SFD. The 
optical depth in the V-band can be found from the colour 
excess E(B-V) maps, 



t(V) 



A{V) 
1.086 



= 2ME(B - V), 



(5) 



where we have us ed a mean galactic value A(V) / E(B — 
V)='6.l (Whittet |l992|) . Reach et al. ( |l995| ) suggest that 
dust emitting in the wavelength range 100-300 /im traces 
interstellar extinction. Since the FIR optical depth in eq. 
(||) has been measured using data at 100 and 240 /im, it 
is then justified to compare it with extinction as in eq. (|5|) 
to find the ratio of the extinction coefficient and emissiv- 
ity. Knowing the optical depths from (||) and (|5J) , we can 
compute a map of the ratio as in eq. (@); we obtain a mean 
value of 



gextQO 

cra (100 /im) 



= 760 ± 60 



for both hemispheres. This value is included, together with 
other multiplicative factors, in the calibration coefficient 
p as in eq. (22) in SFD. As pointed out by the referee, 
an estimate for Q O xt(^0/<3cm can be easily derived from 
that equation, if the DIRBE colour corrections factors, 
slowly depending on T, are omitted. Following this way 
we obtained a value of 765.5. SFD give an error of 8% for 
p and this is the value quoted here. Since most (« 90%) 
of the elliptical galaxies used to calibrate colour excess 
maps have galactic latitude b > 20°, one may argue that 
the measured value is characteristic only of high latitude 
dust. 



Reach et al. ( 1995 ) find that the emissivity (eq. |2j) is 
best determined by fitting the FIRAS spectrum on the 
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Galactic plane. They say that high latitude data have a 
smaller signal-to-noise ratio and can be fitted satisfactory 
with (3 = 2 (eq. [I]) although the same emissivity as on the 
plane cannot be excluded. Under the hypothesis that the 
same kind of dust is responsible for the diffuse emission in 
the whole Galaxy, we have corrected SFD temperatures 
using Reach et al. emissivity (eq. @). The new tempera- 
tures are a few degrees higher than those measured with 
(3 = 2 (as an example we pass from a mean value of 18K 
in a 20° diameter regions around the north pole to a new 
estimate of 21K). It is interesting to note that the differ- 
ence between the two estimates of temperature is of the 
same order as the difference between the temperatures of 
warm dust at high and low Galactic latitude in Reach et 
al. ( 1995 ) and this may only be a result of the different 
emissivity used to retrieve the temperature. 
When the correction is applied 



Qcxt(U) 
Q em (100 Aim) 



= 2390 ±190. 



The new ratio is about three times higher, and this is a 
reflection of the change of temperature in the black body 
emission in (Q): for a higher temperature, a lower emis- 
sivity in the FIR is required to produce the same emis- 
sion. Uncertainties in Q e xt(V)/Q om (\o) are thus greatly 
affected by assumptions about the emissivity spectral be- 
haviour. 



3. Comparison with other measurements 

We now compare our emissivity for (3 = 2 with literature 
results derived under the same hypothesis. Since no emis- 
sivity has been derived to our knowledge assuming eq. (0), 
we do not attempt any comparison with that result. All 
the data are scaled to A = 100 ^m. 

Studying the correlation between gas and dust emis- 
sion from FIRAS and DIRBE, Boulanger et al. ( |l996|) 



derived an emissivity t/Nh = 1.0 ■ 10 cm at 250 
/um for dust at high galactic latitude; assuming the 
canonical Nh = 5.8 • 10 21 E(B — V ) cm -2 mag -1 and 
A(V)/E(B - V)=3.1 (Whittet |1992| ), this is equivalent 
to Q ex t(^)/<9em(100 /im)= 790. Q uite similar values are 
found in the Draine & Lee ( |l984| ) dust model, which has 
a (3 = 2 spectral dependence in this wavelength range. At 
125 /im the optical depth is t/Nh — 4.6 ■ 10~ 25 cm 2 which 
corr espon ds to <3ext(V)/Qcm(100 /zm)=680. Sodroski et 
al. ( 19971 ) finds a value for the ratio at 240 /im, using 
literature data identifying a correlation between B-band 
extinction and 100 /im IRAS surface brightness in high 
latitude clouds, assuming a dust temperature of 18 K. 
Converted to our notation, using a standard extinction 
law, the ratio is Q e xt(^)/Qom(100 /im)=990. 



The measurement by Whitcomb et al. (1981) on the 



reflection nebula NGC 7023 is the mos t com monly quoted 
value for the emissivity (Hildebrand 1985 ). Their value 



8000 - 




1 Mm I 



Fig. 1. FIR emissivity derived in this work: assuming (3=2 
(solid line) and the Reach et al. (1995) spectral depen- 
dence (dotted line). Error bars at 125 ^m give the range 
of values of Whitcomb et al. (1981). The error bar at 250 



fim gives the range of values of Casey (1991). Data points 
are from Draine & Lee (1984) model (st ars), Sodroski et 
al. ( 1997 ) (square) and Boulanger et al. ( 1996 ) (circle). 



derived at 125 /im for (3 = 2 is only marginally consis- 
tent with our result. Following our notation, their result 
is equivalent to Q O xt(U)/Q C m(100 /j,m) = 220 and 8000 
using the energy balance and the extinction meth od, re- 
spectively. The values obtained by Casey ( 1991 ) on a 
sample of five nebulae using the energy balance method 
are a factor of 3 smaller than ours (corresponding to 
Ocxt(U)/Q cm (100 ^m) = 80-400). 

In Fig. |l| we show the literature data (plotted at the 
wavelength they have been derived in the original papers) 
together with our derived emissivity laws. We have added 
the value for Draine & Lee (1984) model at 250 /im. 



4. Gas-to-dust ratio of external spiral galaxies 

We now exploit the FIR emissivity derived in this work by 
determining dust masses for nearby spiral galaxies. Follow- 



ing Hildebrand (1983) dust masses can be measured from 
FIR emission using 



-Wdust = 



4ap 



F{X)D 2 

B(X,T d ) '3Q om (A) 



(6) 



where F(X) is the total flux at the wavelength A, D the 
distance of the object, B(A,Xd) the Planck function, a 
the grain radius (0.1 /im) and p the grain mass density 



1 Whitcomb et al. ( 198l|) and Casey ( 1991 ) originally pre- 
sented values for Q a ^%(UV)/Q CY a{FIR): we have corrected to 
QsKt(V)/Q em (FIR) using the provided t(UV) = 2r(V). 
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Galaxy Gas Mass emissivity (|1|) emissivity (H) 

10 10 M Q T(K) G/D ratio T(K) G/D ratio 



NGC 628 


1.1 


16 (17) 


90 (190) 


18 (20) 


100 (200) 


NGC 660 


0.91 


19 (21) 


110 (230) 


23 (26) 


120 (250) 


NGC 5194 


0.75 


18 (20) 


90 (180) 


21 (24) 


90 (190) 


NGC 5236 


3.5 


19 (21) 


240 (500) 


22 (25) 


255 (540) 


NGC 6946 


3.0 


17 (19) 


75 (150) 


20 (22) 


80 (160) 


NGC 7331 


1.0 


17 (19) 


70 (145) 


20 (22) 


70 (155) 



Table 1. Sample of galaxies from Alton et al. (1998a). Gas masses have been derived from Devereux & Young ( |1990[ ) 
(Van Driel et al. (1995) for NGC 660) and corrected to the distances quoted by Alton ct al. (1998a). Dust 
temperature and gas-to-dust ratio are derived using eq. (|l|) with (3=2 (Q ex t(^)/Qem(100 /Ltm)=760), and eq. (J2|) 
(Qext(V)/Q cm (100 ^m)=2390). Values within brackets are derived under the hypothesis that ISO fluxes are overesti- 
mated by 30%. 



(3 g cm ). The emissivity <5om(A) is derived from t he ra- 

tio Q cxt (V)/Qc 

Whittet]l99f). 



t (A) assuming Qext(V)=1.5 (Casey 1991 



Alton et al. ( 1998a ) provide total fluxes at 100 /im and 
200 fj,m from IRAS and ISO for a sample of spiral galax- 
ies. We have derived dust temperatures and masses using 
Qcxt( V)/Qc m(100 /im)=760 and 2390, for 0=2 and Reach 
et al. ( 1995 ) emissivities, respectively. Using literature val- 
ues for gas masses, we have computed the gas-to-dust ra- 
tios. Values of gas masses, temperatures and gas-to-dust 
ratios are presented in Table |]. 

The mean value of the gas-to dust ratio for the sample 
is 100 using eq. (0), 110 using eq. ph. Mean temperatures 
go from 18K with the = 2 emissivity to 2 IK when the 
Reach et al. ( |1995[ ) behaviour is assumed (as for the north 
galactic pole in Sect. |^). Alton et al. ( |1998a[ ) pointed out 
that ISO 200 fim fluxes could be overestimated by about 
30%; correcting for this we obtain a mean gas-to-dust ratio 
of 220-240 (for (3=2 and Reach et al. Ql995|) emissivity, 



respectively). As shown above, dust masses obtained with 
the two methods are quite similar. This can be explained 
substituting eqs. (Q) and (Q) into (|6|). For A = 100 ^im wc 
can derive 



M d 



5(100 Aim, T d G ) 



5(100 /im,T d ) ' 

where T^ is the mean temperature of dust in the Galaxy. 
From the equation it is clear that the dust mass determi- 
nation is insensitive to the emissivity law used, as long as 
the dust temperature in external galaxies and in our own 
are similar. 

Our range of values for the gas-to-dust ratio (100- 
23 0) enc ompasses the Galactic value of 160 (Sodroski et 
al. 1994[). As a comparison, the mid- value of Whitcomb et 
al. (1981) would have given dust-to-gas ratios larger by a 
factor 1.5. 

5. Conclusion 

We have derived the dust emissivity Q cm in the FIR using 
the wavelength dependence derived from the FIR Galac- 



tic spectrum (Reach et al. 1995 ). The emissivity has been 
normalised to the extinction efficiency in the V band using 
dust column density maps calibrated to Galactic extinc- 
tion (SFD). Q om depends strongly on the assumed wave- 
length dependence. For a (3 = 2 emissivity index we ob- 
tained 



Qcm (A) — 



Qcxt(V ) f 100 /an 
760 



This result is consistent with other values derived from 



FIR Galactic emission (Boulanger et al. 1996; Sodroski et 
al. |1997p and with the Draine & Lee ( |1984| ) dust model . 
The widely quot ed emissivities of Whitcomb et al. (|1981|; 



Hildebrand |1983| ) derived from the reflection nebula NGC 
7023 are only marginally consistent with our values, while 
the emissivity measured by Casey (1991) on a sample of 
five nebulae are smaller by a factor of 3. This may suggest 
a different grain composition for dust in the diffuse inter- 
stellar medium compared to reflection nebulae. 

When the wavelength dependence derived by Reach et 
al. (1995) on the Galactic plane is used, we obtain 



QomW 



<2ext(V) / 100 /an 



2390 



A 



2.005 



1 + (200 nm/Xy 



1/6' 



We have used the derived emissivities to measure dust 
masses from 100 /tin and 200 /im fluxes of a sample of 
six spiral galaxies (Alton et al. 1998a ) . We have retrieved 
similar dust masses with both the spectral dependences. 
The gas-to dust ratios of our sample (100-2 30) a re close 
to the Galactic value of 160 (Sodroski et al. 1994 ). 
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